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Abstract. We review selected results from a recent in-depth study of jet shapes and jet cross sections in 
ultra-relativistic reactions with heavy nuclei at the LHC [1]. We demonstrate that at the highest collider 
energies these observables become feasible as a new, differential and accurate test of the underlying QCD 
theory. Our approach allows for detailed simulations of the experimental acceptance/cuts that help iso- 
late jets emerging from a dense QGP. We show for the first time that the pattern of stimulated gluon 
emission can be correlated with a variable quenching of the jet rates and provide an approximately model- 
independent approach to determining the characteristics of the medium-induced bremsstrahlung spectrum. 
The connection between such cross section attenuation and the in-medium jet shapes is elucidated. 

PACS. 13.87.-a Jets in large-Q^ scattering - 24.85.-|-p Quarks, gluons, and QCD in nuclear reactions - 
12.38.Mh Quark gluon plasma 



1 Introduction 

At present, most measurements of hard processes in heavy 
ion collisions are limited to single particles and particle 
correlations, which are only the leading fragments of a jet. 
There is general agreement that final-state parton energy 
loss in the QGP controls inclusive particle suppression, 
the experimental methodology of determining Raa(j>t) 
or /^yi (pti , PT2 ) is well-established, and the theory of jet 
quenching has been very successful [2j. The data, however, 
cannot resolve the staggering order of magnitude "system- 
atic" discrepancy in the medium properties extracted via 
competing phenomenological models [3]. It is, therefore, 
critical to find alternatives that accurately reflect the en- 
ergy flow in strongly-interacting systems and have a more 
direct connection to the underlying QCD theory. 

The intra-jet energy distribution and the related cross 
section for jets in the case of heavy ion reaction closely 
match the criteria outlined above. In this paper we study 
the magnitude of in-medium modification for these ob- 
servables in Pb-|-Pb collisions at -^/s = 5.5 TeV at LHC. 
We demonstrate that a natural generalization of leading 
particle suppression to jets. 



dyd-'ET 



(1) 



amplifies the observable effect and the jet radius R'^'^^ 
and the minimum particle/tower energy min ~ w™'" 
provide, through the evolution of i?iJ^(£;T; i?'"'''', 
at any fixed centrality, experimental access to the QGP 
response to quark and gluon propagation. 

We refer the reader to [4] for a summary of the com- 
plications in defining a jet and the related topic of jet- 
finding algorithms. To make our discussion simpler, we 
will assume that the complications of the different defini- 
tions can be subsumed into an Rsep parameter [T] . Once a 
jet axis and all of the jet particles / calorimeter towers "i" 
have been identified, the "integral jet shape" is defined as: 



E^iErhOir ~ (R 



■jct)i) 



E,(i^T).e(i?-(i?jct).) 



(2) 



where r, R are Lorentz-invariant opening angles, i?y = 
y^irji — r/j)^ + {(pi — and i represents a sum over all 
the particles in this jet. 'I'int{r;R) is the fraction of the 
total energy of a jet of radius R within a sub-cone of radius 
r. It is automatically normalized so that S'lnt (i?; i?) = 1 
and the differential jet shape is defined as follows: 



i;{r;R) = 



d^int{r]R) 
dr 



(3) 



is sensitive to the nature of the medium-induced energy 
loss. The steepness of the final-state differential spectra 
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Understanding the many-body QCD theory behind jet 
shape calculations will naturally lead to understanding the 
attenuation of jets in reactions with heavy nuclei. 
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2 Theoretical results for vacuum jet shapes 

Jet shapes in vacuum will provide the baseline for jet 
shape studies in nuclear collisions. Here, we follow the 
methodology outlined in Ref. t5j and generalize this ap- 
proach to include finite detector acceptance. 



2.1 Leading order results 

A fixed order QCD calculation of jet shapes is based on 
the splitting functions Pa-,bc{z), the distributions of the 
large fractional lightcone momenta (or approximately the 
energy fractions) of the fragments relative to the parent 
parton, z and 1 — z respectively. To lowest order, recalling 
that ipaif] R) describes the energy flow cx z, we can write: 



a, 2 
2'K r 



dzzPa^bc{z), 
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with Z is defined as follows: 




Z — max ^ Zmim 

Z = max <, ZYnim 



r + R 
r 



if r < {Rsep - l)R , (5) 
if r > {Rsep - l)R . (6) 



In Eq. ^ r — (1 — z)p is related to the opening an- 
gle p between the final-state partons. In "elementary" 
p + p collisions the inclusion of soft particles {zmin ~ 0) 
in theoretical calculations is not a bad approximation. 
In heavy ion reactions, especially for the most interest- 
ing case of central collisions, there is an enormous back- 
ground of soft particles related to the bulk QGP proper- 
ties. Jet studies will likely require minimum particle en- 
ergy > 1 — 2 GeV at RHIC and even more stringent cuts 
at the LHC. Furthermore, control over Zmin can provide 
detailed information about the properties of QGP-induced 
bremsstrahlung. The kinematic constraints on the values 
of z are discussed in detail in [T]. We find for the jet shape 
functions for quarks and gluons: 
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In the Zmin limit Eqs. ([7]) and ([5]) reduce to previously 
known results j5j. 



Fig. 1. Numerical results for the differential jet shape in p-l-p 
collisions at ^/s = 5.5 TeV at the LHC. Solid lines represent 
jet shapes with R = 0.7, ci;""" — GeV, dashed lines stand for 
jet shapes with R = 0.4, lj™'" = GeV, and dashed-dotted 
lines are for jet shapes with R = 0.7, a;™'" = 10 GeV. The 
inserts show integrated jet shapes 'I'int{r; R). 



2.2 The full vacuum jet shape and numerical 
simulations 

A number of important effects beyond the fixed order, 
Eqs. (O and ([8]), contribute to the observed jet shape [1] 
[S]. Taking all contributions and ensuring that there is no 
double counting at small r/R to 0{as) we find: 

V-W = V-coll W (Pir) - 1) + ^LD W + V-^.LD W 

+V^pc(r)-f V.,PCW : (9) 

On the right-hand-side of Eq. ([9]) the first term comes 
from Sudakov resummation with subtraction of the lead- 
ing 1/r, (1/r) log(l/r) contribution at small r/R to avoid 
double counting with the fixed order component of the 
differential jet shape. The second and third terms repre- 
sent the leading-order contributions in the final-state and 
the initial-state. The last two terms represent the effect of 
power corrections. In a full calculation the relative quark 
and gluon fractions fq+fg = 1 are also needed: ^p{r; Et) = 
fg{ET, ET)+fg{ET, \/i)^g(r; Et). This theoret- 

ical model is validated against CDF II data [1] . 

Figure [T] shows our numerical results for the jet shape 
for two different energies Et = 50 GeV and 250 GeV 
and two cone radii R — 0.7, 0.4 in p + p colhsions at 
^/s = 5.5 TeV at LHC. When plotted against the rela- 
tive opening angle r/R, these shapes are self-similar, i.e. 
approximately independent of the absolute cone radius 
R. One of the main theoretical developments in Ref [1] 
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is the analytic approach to studying finite detector ac- 
ceptance effects or experimentally imposed low momen- 
tum cuts. In Fig. [T] this is illustrated via the selection of 

Zmin = PT min/ Et = 0.2, 0.04 {PT min = w"'" = 10 GcV) . 

Eliminating the soft partons naturally leads to a narrower 
branching pattern. However, for this effect to be readily 
observable 10 — 20% of the jet energy, going into soft par- 
ticles, must be missed. Thus, even with j>x' min 

-few GeV 

cuts in Pb-|-Pb collisions at the LHC aimed at reducing or 
eliminating the background of bulk QGP particles that ac- 
cidentally fall within the jet cone, the alteration oiip{r/R) 
is expected to be small. The integral jet shape 'Pintir/R), 
shown in the inserts of Fig. [U is not optimal as a tool for 
identifying kinematic and dynamic effects on jets due to 
its much smaller sensitivity. 



Jet axis 




6)™"(2) 



3 Medium-induced contribution to the jet 
shape 

The principal medium-induced contribution to a jet shape 
comes from the radiation pattern of fast quarks or glu- 
ons, stimulated by their propagation and interaction in the 
QGP. There is a simple heuristic argument which allows 
one to understand how interference and coherence effects 
in QCD amplify the difference between the energy distri- 
bution in a vacuum jet and the in-medium jet shape [6]. 
Any destructive effect on the integral average parton en- 
ergy loss AE^"''^, such as the Landau-Pomeranchuk-Migdal 
effect, can be traced at a differential level to the attenu- 
ation or full cancellation of the coUinear, fc^ <C l^j, gluon 
bremsstrahlung : 



A prad 

^-'^LPM suppressed 



dI3 



LPM suppressed 



dujd^kT 



[kr < w)lpm suppressed 



and we indicate the parts of phase space where the modi- 
fication of the incoherent dl^ / duid^kx is most effective. 

In our calculation we use the GLV formalism of ex- 
panding the medium-induced radiation in the correlations 
between multiple scattering centers [71[5]. Theoretical in- 
terest in the angular bremsstrahlung distribution was stim- 
ulated by experimental measurements of enhanced trig- 
gered opposite-side particle correlations away from the 
naive Z\(/) = tt pifTUlfTT] . It has been know that the leading 
n = 1 contribution to final-state medium-induced radia- 
tion does not have a collinear with the jet axis compo- 
nent [12]. We now show [I] that this result is general and 
holds to any order in the opacity expansion: 



hm ^ 



1-^0 duidcfidr 



(11) 



Numerical simulations, using Monte-Carlo techniques, con- 
firm independently that dl^ / diod^'k vanishes as k 

. . . 

The implication of our finding is that there is very lit- 
tle overlap between the techniques used to compute the 



Fig. 2. Schematic illustration of the cone radius R and the 
particle/tower pt / Et selection. The measured energy is the 
one that comes from particles with pt > uJmin and within R. 

"vacuum" and medium-induced contributions to the jet 
shape. Also, selecting different jet radii R and pt min of 
the particles will significantly alter both the jet shape and 
the amount of energy lost by the hard parton which can 
be recovered in the experimental measurement. A clear 
strategy will be to use the leverage arms provided by 
R Ji"^°-^ in the evaluation of the Z\i?rad) and pT min 
(= u;™™ in the evaluation of the AEn^^) to determine the 
distribution of the lost energy. This is illustrated schemat- 
ically in Fig. [51 Theoretically, the first quantity to be cal- 
culated is: 



(10) AE'' 
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The separate dependence of AE^'^{R 
the cone radius and the momentum acceptance cut is 
more clearly illustrated in Fig. [3l We show central, mid- 
central and peripheral collisions, impact parameters b = 
3, 8, 13 fm, respectively, in Ph 4- Pb reactions at LHC 
at nominal y^s. We notice that, not surprisingly, the ratio 
AE™{R™'^^, 10"^™) /E goes down at larger impact parame- 
ters because the energy loss of the jet decreases in periph- 
eral collisions. More importantly, at each impact parame- 
ter there is a variation of the amount of the bremsstrahlung 
energy, recovered in the cone. This is precisely the varia- 



tion that will map on the R\j^ {Et, 
able. 



^) observ- 



4 Tomography of jets in heavy ion collisions 

The purpose of this Section is to relate the theory of jet 
propagation in the QGP to experimentally measurable 
quantities. 
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Fig. 3. 2D projections of Eq. (fT2|l . The left panel shows 
the fractional energy loss dependence on the jet radius 
(a;°"° = 0) and the right panel shows this dependence ver- 
sus the acceptance cut ti^""" {'R^'^^ = A gluon jet of 
Ejet = 20 GeV in 6 = 3, 8, 13 fm Pb + Pb collisions at LHC 
was used as an example. 



4.1 Experimental observables 

An essential ingredient that controls the relative contri- 
bution of ipva.c.{r/R) and ipmedXr/R) = {1/ AE,^^)dIs / dr 
to the observed differential jet shape in heavy ion reac- 
tions and also determines the attenuation of the jet cross 
sections is: 



/ 



^ zAii;,ad{(o,fl);(^"'",i?)} 

Z\£;,ad{(0,i?-);(0,S)} 



(13) 



the fraction of the lost energy that falls within the jet 
cone, r < R, and carried by gluons oiuj > w'"'" relative to 
the total parton energy loss without the above kinematic 
constraints. If this fraction is known together with the 
probability distribution P(e) for the parton energy loss 
the medium-modified jet cross section per binary N + N 
scattering can be calculated as follows: 
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NN 
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d^E'dy 



(14) 



where E'.^ = Et/{1 - (1 - /,,g) • e). The (1 - /,,g) • e fac- 
tor accounts for the total "missed" energy in a jet cone 
measurement, which necessitates Elp > Et- Next, we ob- 
tain the full jet shape, including the contributions from 
the vacuum and the medium-induced bremsstrahlung: 
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_ET-dependent nuclear 



modification factor 
(^max^^mm^ different jet cone radii 7?'"=''' at 

6 = 3 fm in central Pb + Pb collisions with ^/s = 5500 GeV. 



normalized to unity. Integrating over r in Eq. (|15p , we can 
easily see that the correct "Norm" is the quenched cross 
section, Eq. (fT4| . Discussion of simple limiting cases for 
Eqs. (fTi|) and which provide insight into the contri- 
butions to the full jet shape and the jet cross section is 
given in [1]. 



4.2 Numerical results 



Clearly, the centrality dependence of R-^^ (i?™'''', w™'") 
at different impact parameters will remain an important 
handle in understanding the QCD medium. But the ma- 
jor development, that we report on, is the emerging un- 
derstanding on how to get directly at the mechanisms of 
jet interaction in matter. In Fig. [5] we show the evolution 
of the nuclear modification factor for the jet cross section 
by varying the size of the cone radius R"^^^ (w™'"^'^). We 
see that with increasing cone radius R™^^ the quenching 
of the jet cross section disappears at all Et and finally 
reaches unity when R'^'^^ = 2.0 and all of the lost en- 
ergy is re-captured. It is important to note the factor of 5 
to 10 variation in the quenching of the jet cross section. 
Such variation will pinpoint the angular distribution of 
the bremsstrahlung gluons. We have also studied in detail 

the dependence of i?^^ (i?™'''^, w™™) on the acceptance 
cut w™™ [T] and found similar sensitivity to the gluon en- 
ergy dependence of dl^/dio. The continuous variation of 
quenching values may help differentiate between compet- 
ing models of parton energy loss, thereby eliminating the 
order of magnitude uncertainty in the extraction of the 
QGP density 

Next, we turn to the pattern of energy flow for in- 
medium jets. In Fig.Othe jet shape in vacuum, the medium- 
induced jet shape and the total jet shape in the QGP 
with different selection of cone radii and different trans- 
verse energies Et in central Pb + Pb collisions at y/s = 
5.5 TeV are shown together for comparison. An interest- 
ing conclusion is that there is no significant distinction 



[(1 - 6) rjc. (r/R) 

+ Ua-^^^r^d.(r/R)] .(15) 



We recall that, by definition, the area under any differen- 
tial jet shape, -0101. (r/R), V'vac. (r/R) and TAmod. (r/R), is 
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Fig. 5. Comparisons of the jet shape in vacuum, medium- 
induced jet shape, and the total jet shape for cone radii 
R = 0.7, 0.4 and for different jet energies Et — 
20, 50, 100, 200 GeV, respectively, in Pb + Pb collisions at LHC. 



Table 1. Summary of the mean relative jet radii in vacuum 
and in the QGP medium. Jet cone radii R = 0.4 and R — 0.7 
and energies Et ~ 20 and 100 GeV at -^s = 5.5 TeV central 
Pb + Pb collisions at the LHC are shown. 
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Fig. 6. The ratios of the total jet shape in heavy-ion col- 
lisions to the jet shape in vacuum for jet energies Et ~ 
20, 50, 100, 200 GeV and cone radius 7? = 0.4 in fe = 3 fm 
Pb + Pb collision at = 5.5 TeV. 



Finally, we point out where the effect of the QGP 
medium can be more readily observed in the differen- 
tial jet shapes. These are the "tails" of the energy flow 
distribution in tptotXi" / R), near the periphery of the jet 
cone. We can see that, with cone radius R = 0.4, the 
ratio tptot.{^ / R) / '>pvac.{r / R) could reach about 1.75 when 
r/R — > 1. However, for experiments to reliably detect and 
quantify this enhancement of the ratio of the total jet 
shape in the QCD medium to the jet shape in the vacuum 
tn the r/R > 0.5 region, high statistics measurements will 
be necessary. The required precision can be achieved with 
the large acceptance experiments at the LHC beyond the 
proof-of-principle jet identification at RHIC [H]. 



between the jet shape in vacuum and the total in-medium 
ipir/R). The underlining reason for this surprising result is 
that although medium-induced gluon radiation produces 
a broader V'mod.(?'/^), this effect is offset by the fact that 
the jets lose a finite amount of their energy. Furthermore, 
when part of the lost energy is missed due to finite accep- 
tance, the required higher initial virtuality jets are inher- 
ently narrower. 

In Table[T]we present the mean relative jet radii {r/R) 
in the vacuum and in the QGP medium created at the 
LHC for two different cone selections R — 0.4 and R = 0.7 
and two jet energies Et = 20, 100 GeV. We see that in the 
realistic numerical simulation of quark and gluon propa- 
gation in the QGP there is very little < -f 10% variation 
in this observable. This characteristic difference is slightly 
larger for a smaller cone, since it emphasizes the large- 
angle character of the medium- induced radiation [7lfT2] . It 
is important to stress that the QGP is rather "gray" than 
"black" and only a fraction of the energy of the leading 
jet is lost in the medium. The effect of even a moderate 
AE^^^max^ ^min-j ^^^^ ampfificd by the steeply faUing 

cross sections for the R'^^ {Et', R™^^,i^™^'^) observable 
but this is not the case for {r/R). 



5 Conclusions 

The unprecedentedly high center of mass energy at the 
LHC will open a new frontier for jet physics in dense QCD 
matter: tomography o/jets that can help pinpoint energy 
ffow in strongly-interacting systems and significantly re- 
duce uncertainties in the QGP properties extracted from 
different phenomenological models. In this summary we 
highlighted several important results from a recent com- 
prehensive study of jet shapes and jet cross sections in 
y/s = 5.5 TeV Pb + Pb collisions at the LHC JLj. A the- 
ory of jet shapes in the vacuum, generalized to deal with 
experimental acceptance cuts, was discussed and then ap- 
plied to study jet shapes in heavy-ion collisions by in- 
cluding the contribution from the medium-induced gluon 
bremsstrahlung. Realistic numerical simulations of jet cross 
sections and jet shapes were used to illustrate new sig- 
natures of quark and gluon propagation in the QGP. It 
was shown that a suitable generalization of the nuclear 

modification factor for jets, Rj^j^ {R'^^^^, lj™™), provides a 
sensitive tool to probe the energy and angular dependence 
of the stimulated gluon emission pattern. The largest dif- 
ference between the jet shape in vacuum and the total 
in-medium jet shape was manifest in the periphery of the 
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cone, r/R 1, and for smaller radii, e.g. — 0.4. The 

enhancement of the mean relative jet radius {r/R) due to 
QGP effects was rather modest. 
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